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We develop a nonstandard concept of atomic clocks where the blackbody radiation shift (BBRS) 
and its temperature fluctuations can be dramatically suppressed (by one to three orders of mag- 
nitude) independent of the environmental temperature. The suppression is based on the fact that 
in a system with two accessible clock transitions (with frequencies V\ and u 2 ) which are exposed 
to the same thermal environment, there exists a "synthetic" frequency u Byn tx {y\ — £12^2) largely 
immune to the BBRS. As an example, it is shown that in the case of 171 Yb + it is possible to create 
a clock in which the BBRS can be suppressed to the fractional level of 10" 18 in a broad interval 
near room temperature (300±15 K). We also propose a realization of our method with the use of 
an optical frequency comb generator stabilized to both frequencies v 1 and v 2 - Here the frequency 
K; yn is generated as one of the components of the comb spectrum and can be used as an atomic 
standard. 

PACS numbers: 03.75.Dg, 06.20.F-, 37.25. +k, 42.62.Fi 

The main progress in modern fundamental metrology 
is connected with the development of atomic clocks. The 
most promising frequency standards today are based on 
single trapped ions [l( and on ensembles of neutral atoms 
confined to an optical lattice at the magic wavelength 
[H, HJ. It is believed that these clocks can provide fre- 
quency references with unprecedented small systematic 
uncertainties in the 10~ 17 -10~ 18 range. This progress 
will probably lead to a redefinition of the unit of Time 
and to new fundamental tests of physical theories in par- 
ticular in the fields of General Relativity, cosmology, and 
unification of the fundamental interactions 0, [1| . 

The largest effect that contributes to the systematic 
uncertainty of many atomic clocks is the interaction of 
the thermal blackbody radiation with the atomic eigen- 
states. This effect was first considered in 1982 for cesium 
atomic clocks but remains up to now a major prob- 
lem for many modern atomic time and frequency stan- 
dards. At present there exist three approaches to tackle 
the blackbody radiation shift (BBRS) problem. The first 
one is the use of cryogenic techniques to suppress this 
shift to a negligible level. This approach is pursued for 
the mercury ion clock Q, for the Cs fountain clock Q, 
and for the Sr optical lattice clock [§]. The second ap- 
proach is the precise temperature stabilization of the ex- 
perimental setup in combination with theoretical and/or 
semicmpirical numerical calculations of the shift at given 
temperature [l(J 111- The third approach is based on the 
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choice of an atom or ion where both levels of the refer- 
ence transition have approximately the same BBRS. Here 
the most promising candidate is 27 A1 + with a fractional 
BBRS of the reference transition frequency of ~10~ 17 
[E 03, followed by 115 In+ [H, El- However, the lat- 
ter approach limits the choice of candidates for tests of 
fundamental theories. 

In the present paper we propose an alternative method 
allowing us to suppress the BBRS and its fluctuations 
in atomic frequency standards by one to three orders of 
magnitude without using cryogenic techniques and pre- 
cise temperature stabilization. Our approach is based 
on the use of two reference transitions in an identical 
thermal environment. We show that in such a system 
there exists a combined frequency for which the BBRS is 
significantly suppressed over a wide temperature range. 
For instance, a trapped 171 Yb + ion meets this condi- 
tion in a straightforward way, because 171 Yb + has at 
least three suitable reference transitions: an electric- 
quadrupole and an electric-octupole optical transition 
3-171, and a magnetic-dipole radiofrequency (rf) transi- 
tion between the ground-state hyperfine sublevels. Apart 
from laboratory standards, the proposed method can be 
particularly useful in cases where it is impossible to con- 
trol the environmental temperature with sufficient accu- 
racy or to use cryogenic techniques, for instance in trans- 
portable frequency standards or in space-based clocks 
that approach the Sun in order to test the local position 
invariance underlying General Relativity Q ■ 

Our approach is based on the fact that for the large ma- 
jority of transitions in atoms or ions that are of interest 
as frequency standard reference transitions, the temper- 
ature dependence A(T) of the BBRS is very well approx- 
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imated by the law oc T 4 . Consider now two clock transi- 
tions with frequencies and v 2 exposed to the same 
thermal environment, i.e., located in the same probe vol- 
ume. We assume that z^ ' < v 2 °^ . The effect of the 
BBRS on each transition frequency can be represented 
as: 
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where a,j is an individual characteristic of the transition 
j determined by the atomic structure and T is the mean 
temperature of the clock operation. Let us introduce the 
coefficient e 12 = a l /a 2 . As is easily seen, the following 
superposition does not experience the BBRS: VyiT) — 
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In compliance with this we 



define a new "synthetic" frequency v sya as 



h2 V 2 (T)] = R[v\ 
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where R is some numerical multiplier whose value can 
be chosen freely. Thus, one can use the frequency z^ syn 
as a new clock output frequency which is immune to the 
BBRS and to fluctuations in the operating temperature, 
while the thermal shifts c^T 4 of the working frequencies 
i/j can be large. 

One possibility is to independently measure both fre- 
quencies v Y 2 (T) and to use for the clock operation the 
synthetic frequency according to Eq. ((2]) (assuming, for 
example, R = ±1 ). In this case, obviously the synthetic 
frequency does not directly correspond to any frequency 
of a physical signal. Another approach is to synthe- 
size this frequency as a real physical signal by means 
of an optical frequency comb generator. Let us consider 
the situation where two modes of the frequency comb 
generator are stabilized to the two optical frequencies 
v \ 2(^0 = /o + % 2/r at a given temperature T (sec 
FiglTJ) . As a result, the parameters of the comb spectrum, 
i.e., the pulse repetition rate f r and the offset frequency 
fo are unambiguously determined and the frequency of 
the m-th mode equals: 
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From this expression one can define a number m = m 
for which the coefficient of the temperature-dependent 
term is zero: 
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FIG. 1: Illustration of femtosecond comb stabilized to the 
two clock transitions with frequencies Vi and v 2 at a given 
temperature T. 

quency v mo is the synthetic frequency defined in Eq. ^ , 
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as it should be. Here, mo is the natural number closest 
to the value of the right-hand side of Eq. For this 
it is necessary to satisfy the condition m > that is 
equivalent to ^iy™*' > in Eq. (J5J). 

Apart from the frequency v ma which is a component 
of the optical spectrum of the frequency comb genera- 
tor, in our system one can also define the much smaller 
frequency 



fr + 
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which corresponds to a rf standard at f mo /m . Since the 
frequencies f r and f can be extracted from a stabilized 
comb generator with negligible error, one can use them 
to synthesize v mo /m . This synthesized radiofrequency 
has the same immunity to BBRS as v rno ■ It is interesting 
to note that the radiofrequency given in Eq. ((5]) is well- 
defined in our system even if m < in Eq. ([4]), i.e., if 
the basic frequency component v mo exists only virtually. 

As was shown above, in the case of a frequency comb 
stabilized to two BBR-shiftcd clock transitions with fre- 
quencies V\(T) and v 2 (T), there exists a frequency com- 
ponent v mo (for to > 0) for which the thermal shift and 
the sensitivity to temperature fluctuations vanish. This 
frequency component can serve as an atomic frequency 
standard. In practice, the BBRS is strongly suppressed 
for a range of comb frequencies around v mo . The residual 
shift of frequency components v mo ±i near v ma is given by: 



mo±l 
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This indicates that the suppression is effective as long 
This shows that the BBRS is suppressed for the frequency as (n 2 - n x ) > I. For example, for frequencies vf> and 
v ma . After a simple transformation we see that the fre- v 2 ^ in the optical range, the comb mode index difference 



3 



(n 2 — n-J will typically be in the range of 10 5 or higher 
(see the discussion for the case of 171 Yb + below). On 
the whole, the choice of the exact value of the synthetic 
frequency i/ syn is to some extent arbitrary if one takes 
into account that the coefficient e 12 is only known with 
limited accuracy and that Eq. (1) is an approximation 
that neglects higher-order terms in the temperature de- 
pendence of the BBRS [HI . Including higher-order terms 
the BBRS can be expresssed as: 



F=2 
F=\ 



AW(T) 



+ ... (j - 1,2). (8) 



From this we can estimate a basic limitation of the pos- 
sibility to suppress the BBRS and its temperature de- 
pendence. Usually, near room temperature To = 300 K 
the contribution of the higher terms [6^(T/T ) 6 + •■■] is a 
factor of 10 to 10 3 smaller than that of the main T -term 
@, Ell • This indicates that here it would not be useful to 
suppress the T 4 -dependence of v syn to better than one to 
three orders of magnitude because higher-order contribu- 
tions to the BBRS remain uncompensated. For example, 
in order to achieve a suppression factor of 10 2 for the In- 
dependence, it would be sufficient to know the coefficient 
e 12 with relative uncertainty of 10 -2 . 

It may be noted that apart from theoretical calcula- 
tions the coefficient e 12 can be determined by purely ex- 
perimental means. To do this we can apply a quasistatic 
electric field (or the field of an infrared laser) to determine 
the shifts of the reference transition frequencies v 1 and 
v 2 due to the differences in the static polarizabilities of 
the involved atomic energy levels. From a practical point 
of view it is very advantageous that we do not need to 
know the magnitude of the electric field at the place of 
the atoms because we only have to determine the ratio 
a±/a2- If a frequency comb generator is stabilized to v 1 
and v 2 as shown in Fig.l, the frequency v mo can be iden- 
tified in a direct way as the frequency component which 
does not experience any scalar Stark shift in the applied 
quasistatic field. 

As an example that permits the practical realization of 
the ideas presented above, we consider the ion 171 Yb + . 
As shown in Figj2j the level system of 171 Yb + pro- 
vides two narrow-linewidth transitions from the ground 
state in the visible spectral range which can be used as 
reference transitions of an optical frequency standard: 
the quadrupole transition 2 Si/ 2 (F = 0) -^ 2 D 3 / 2 (T = 
2), A w 436 nm and the octupole transition 2 Si/ 2 (T = 
0) ^ 2 F 7/2 (F = 3), A w 467 nm. More detailed infor- 
mation on the spectroscopy of these transitions can be 
found in p~5l - f]~7L l2fj| . It may be noted that the case of 
171 Yb + is especially attractive because here both clock 
transitions lie in a technically convenient frequency range 
and experience exactly the same thermal environment if 
probed in one ion. 

The BBRS of the quadrupole and octupole tran- 
sitions of Yb + were calculated in Rcf. 2l|. This 
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FIG. 2: (Color online) Section of the energy level scheme of 
171 Yb + , showing the hyperfine levels of the 2 Si/2 ground state 
and the two lowest-lying excited states, which are metastable. 
Hyperfine splittings are not drawn to scale. 



calculation is based on calculated oscillator strengths 
and experimental lifetime and polarizability data. The 
room-temperature BBRS of the quadrupole transition 
is calculated as a qua( j ~ —0.35(7) Hz (fractional shift 
5.1(1.1) x 10 -16 ) and that of the octupole transition as 
a oct = -0.15(7)) Hz (fractional shift 2.4(1.1) x 10" 16 ). 
The relatively small value of a oc t and the large relative 
uncertainty is due to nearly equal shifts of the Si/ 2 and 
2 F 7 / 2 levels. 

Using the results of Ref. [U, for 171 Yb+ we find that 
£ u = a oct/aquad = 0.43(22). We expect that the large 
uncertainty of this value can be reduced to less than 1% 
by improved atomic structure calculations or by a direct 
measurement of e 12 as discussed above. In the following, 
we will not take into account the present uncertainty be- 
cause our conclusions remain qualitatively unchanged for 
all values of e 12 in this uncertainty range. In particular, 
we find the synthetic frequency f S y£ mb ^ ~ 607 THz, corre- 
sponding to a wavelength A syn « 494 nm. This frequency 
lies sufficiently close to the initial reference transitions at 
436 nm and 467 nm that it can be generated as a spec- 
tral component of a femtosecond comb generator that is 
locked to the reference transitions as shown in Fig. 1. 

The higher-order contributions in Eq. $5$i to the BBRS 
of the octupole reference transition are negligible com- 
pared to that of the quadrupole transition. For the lat- 
ter, we find b/a w 0.1 at To = 300 K. As a result, we 
estimate that the BBRS can be suppressed to the frac- 
tional level of 2.7 x 10~ 17 at 300 K with variations at the 
level of ±5 x 10 -18 if the ambient temperature varies in 
a broad interval of ±15 K. 

It is also possible to estimate the frequency interval 



round 



(comb) 



where the components of the comb spec- 



trum have a similar level of suppression of the thermal 
shift and of its fluctuations. For a suppression factor of 
10 2 , this interval has a width of the order of 1000 GHz. 
Thus, for d ~ 100 MHz, the indicated interval contains 
10 4 comb modes, each of which could be used as a stable 
frequency reference. 
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Other variants of BBRS-frce optical frequency stan- 
dards at a synthethic frequency can be conceived based 
on transitions 1 So -^ 3 Po in alkaline-earth-like neutral 
atoms confined in an optical lattice. Consider, for in- 
stance, the combination of the reference transitions of 
strontium (v 1 ps 429 THz, A w 698 ran) and ytter- 
bium (y 2 ~ 518 THz, A w 578 ran). Using the calcu- 



lations in Ref. 



in this case we obtain e 



and an estimated synthetic frequency v. 



(comb) 
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1.69 
648 THz 

(A syn ~ 463 nm). The technical realization of this vari- 
ant requires the operation of two lattice-based clocks with 
different atoms (Sr and Yb) in the same vacuum cham- 
ber. 

So far we have considered examples where both refer- 
ence transition frequencies lie in the optical region. How- 
ever, it is also possible to realize schemes where the high 
frequency v 2 is optical, but the lower frequency v 1 corre- 
sponds to a fine- or hyperfine-structure splitting so that 
it lies in the terahertz or microwave range. In contrast 
to the above example of ion Yb + , which seems unique 
because it provides two optical reference transitions, in 
this case one can find many appropriate schemes that use 
a single atomic species. Also the ion 171 Yb + offers the 
possibility of using the ground-state hyperfinc transition 
F = -> F = 1 at ^=12.6 GHz (see Fig© as a low- 
frequency reference transition. For the combination with 
the octupole transition 2 Si/ 2 (F = 0) — > 2 F 7 / 2 (F = 3) at 
1^2=642 THz, a numerical estimate yields e 12 «6.6xl0 -5 
and a synthetic frequency ^ syn = — (v l — £12^2) ~30 GHz. 
(We expect that our estimate on e 12 is accurate to ±20 %, 
but it is not principal for further results.) Here the 
T 6 -contribution to the BBRS (sec Eq.©) limits the 
BBRS suppression at the fractional level of 7.5xl0~ 19 
at T=300 K with variations of ±2xl0 -19 in the temper- 
ature interval of 300±15 K. However, we should also take 
into account the shift of the ground-state hyperfine lev- 
els (oc T 2 ) due to the magnetic blackbody radiation field 
Q . The corresponding BBRS of the hyperfine frequency 
F = -> F = 1 for m Yb+ is: 



A« (T) = -1.616 x 10 



\ 300 / 



Hz. 



(9) 



For f syn =30 GHz this shift results in a fractional level of 
5.4X10" 18 at T=300 K with a variation of ±5xl0" 19 for 
300±15 K. Since the magnetic BBRS contribution can 
be readily calculated with an accuracy of less than 1%, 
it is possible to apply a corresponding correction to z/ syn 
with an uncertainty contribution of less than 10 -19 . 

For 171 Yb + we have thus shown the possibility to 
create a synthetic-frequency-based atomic clock with 
a fractional uncertainty contribution due to BBRS of 
<1.5xl0~ 18 in a broad interval of 300±15 K. To achieve 
this, we need to know the coefficient e 12 with a relative 
accuracy in the range of 0.1-0.2%. In order to reduce the 
BBRS uncertainty contribution to less than 10 -17 , the 



value of e 12 needs only be known with a relative accu- 
racy of 3%. We also have pointed out that for 171 Yb + the 
combination of the octupole optical clock transition with 
the ground-state hyperfine transition can yield a much 
better BBRS suppression than the combination of the 
octupole and quadrupole optical clock transitions. The 
use of the quadrupole transition yields a lower BBRS 
suppression because the upper level 2 D 3 / 2 is connected to 
the 2 Pi/2 level by a strong infrared transition at 2.44 /xm, 
which produces a relatively large T 6 -contribution to the 
BBRS. The final comparison of the two options for BBRS 
suppression should also include detailed estimates on the 
magnitudes of other systematic uncertainty contributions 
in the considered experimental setup. 

The concept of a synthetic atomic frequency standard 
based on two reference transitions can also be extended 
to the case that both reference frequencies isi. 2 lie in 
the microwave range. Atomic fountain clocks are based 
on reference transitions in the microwave range between 
the ground-state hyperfine sublevels of alkali atoms. For 
a synthetic atomic fountain frequency standard, for in- 
stance the combination 87 Rb [y x « 6.8 GHz) and 133 Cs 
(y 2 w 9.2 GHz) can be considered. Here, at the synthetic 
frequency (v l — £12^2) ~ 1-9 GHz it is possible to sup- 
press the fractional BBRS of the individual standards by 
two orders of magnitude. It is interesting to note that 
nearly optimal conditions for the efficient suppression of 
the BBRS are realized in the dual Rb/Cs fountain clock 
described in Ref. [23| because here both reference tran- 
sitions are exposed to the same thermal environment. 

We finally note that the 171 Yb + optical frequency stan- 
dard is a very sensitive system for a search for tempo- 
ral variations of the fine structure constant a 
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The frequencies of the quadrupole and octupole refer- 
ence transitions of Yb + have significant contributions 
from relativistic effects and would undergo changes with 
different sign in consequence of a change of a. The syn- 
thetic frequency that eliminates the BBRS retains this 
sensitivity. The a-dependence of an atomic transition 
frequency may be expressed generally as v = -\- qx, 
where x = (a/ao) 2 — 1, vq defines the frequency at the 
present value of the fine structure constant, ao, and q 
quantifies the sensitivity to changes of a [24|. The q pa- 
rameter for the synthetic frequency is simply given by 
9syn = {qi - £10(72)7(1 — Ei,). For Yb+, with q parame- 



ters as given in Ref. [24(, q syn amounts to about -3220 
THz. Comparison with the Yb + synthetic frequency 
I'syn" 1 ^ 607 THz indicates the strong sensitivity. In 
a test for variations of a, the synthetic frequency would 
have to be compared to an "anchor" reference transition 
with small q value, like the 1 So — > 3 Po transition in Al + . 

In conclusion, we have proposed and developed the 
concept of an atomic frequency standard where the fre- 
quency shift due to the ambient blackbody radiation and 
related fluctuations of the output frequency can be sup- 
pressed by one to three orders of magnitude without us- 
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ing cryogenic techniques. We also expect that our re- 
sults will stimulate refined atomic structure calculations 
on Yb + and other atomic systems that are of interest in 
this context. Such calculations can yield precise values 
for the frequency synthesis parameter e 12 and determine 
limitations of the achievable BBRS suppression. 
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